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ABSTRACT 

The emission mechanisms of the blazar 3C 279 are studied by solving the 
kinetic equations of electrons and photons in a relativistically moving blob. The 
7-ray spectral energy distribution (SED) is fitted by inverse Compton scattering 
of external photons. The bulk Lorentz factor of the emitting blob is found to 
be ~ 25, and the magnetic field is found to be ~ 0.3 G. GeV 7-rays are well 
explained by inefficiently cooled electrons because of the Klein-Nishina effects. 
The electron spectrum is not a broken power law with a steeper spectrum above 
a break energy, which is often used to fit the observed SED. The kinetic energy 
density of the nonthermal electrons dominates the magnetic energy density; this 
result is qualitatively the same as that for TeV blazars such as Mrk 421 and Mrk 
501. The 7-ray luminosity of 3C 279 is often observed to increase rapidly. We 
show that one of the better sampled 7-ray flares can be explained by the internal 
shock model. 

Subject headings: quasars: general — quasars: individual (3C 279) — gamma 
rays: theory - radiation mechanisms: nonthermal 

1. Introduction 

Blazars are active galactic nuclei characterized by high energy emission and short time 
variability (e.g., Krolik 1999). Recent observations of blazars show that they are powerful 
sources of high energy emission (Ulrich, Maraschi, & Urry 1997, for review), which is 
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explained by nonthermal emission by particles in relativistic jets. The emission of blazars 
in the GeV energy range was observed by EGRET onboard the Compton Gamma Ray 
Observatory (Thompson et al. 1995; Mukherjee et al. 1997; Hartman et al. 2001). These 
observations show that a large fraction of emission power from blazars is in 7-ray band (see 
Ghisellini et al. 1998 and Collmar 2001 for reviews). 

Multiwavelength observations of 3C 279 (redshift z = 0.538) were performed several 
times (Maraschi et al. 1994; Wehrle et al. 1998; Hartman et al. 2001). Based on these 
results, 7-rays from 3C 279 are found to be explained by inverse Compton scattering of 
external UV radiation (e.g., Inoue & Takahara 1996). The sources of seed soft photons 
(e.g., accretion disks (Dermer & Schlickeiser 1993), broad-line regions (Sikora, Begelman, 
& Rees 1994), and beamed radiation reflected by broad-line region clouds (Ghisellini & 
Madau 1996)) have been proposed by many authors. Recently, Hartman et al. (2001) 
assumed soft photons from both the accretion disk and the broad-line region in order to 
model the emission spectra from 3C 279. Ballo et al. (2002), on the other hand, applied a 
broad-line region model to 3C 279. 

Most emission models have assumed that electrons obey a power law (e.g., Hartman et 
al. 2001) or a broken power law (e.g., Inoue & Takahara 1996; Ballo et al. 2002). However, 
an alternative electron spectrum is possible as a result of the Klein-Nishina (K-N) effects; the 
smaller cooling rate due to the K-N effects makes the electron spectrum harder (Dermer & 
Atoyan 2002), and the emissivity of scattered photons is smaller (Georganopoulos, Kirk, & 
Mastichiadis 2001). Thus, the electron spectrum must be solved self-consistently, including 
the cooling rate in the K-N regime. 

In previous work, we solved the kinetic equations of electrons and photons simultane- 
ously in a relativistically moving blob (Li & Kusunose 2000; Kino, Takahara, & Kusunose 
2002, hereafter KTK). Here we include the inverse Compton scattering of external soft pho- 
tons (external Compton scattering, hereafter EC) and obtain the physical quantities for 3C 
279 such as the bulk Lorentz factor, the energy densities of electrons and magnetic fields, 
etc. 

The variability of blazars on small timescales is often observed (e.g., Wehrle et al. 1998). 
The time variation of fluxes is not only useful for obtaining constraints on the size of the 
emission region but is also important for studying the emission and particle acceleration 
mechanisms (Mastichiadis & Kirk 1997; Kusunose, Takahara, & Li 2000; Li & Kusunose 
2000; Bottcher and Chiang 2002). Flares observed in X-ray and 7-ray regimes might be 
explained by the internal collisions of plasma blobs in a jet and the subsequent cooling. 
Without considering the details of particle acceleration after the collision, it is possible to 
obtain the emission spectrum by assuming that the time evolution occurs quasi-steadily. In 
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this Letter we present a model to explain the flare of 3C 279 observed in 1996 February. 

In §2, kinetic equations for electrons and photons are described. Numerical results are 
given in §3 for a steady state. Results for time variability are shown in §4. Finally, §5 is 
devoted to the summary of our results and discussion. 



2. Model 

The formulation used in this Letter is the same as in Li & Kusunose (2000) and KTK. 
We assume that a spherical blob with radius R moves with a relativistic speed with Lorentz 
factor T. Accelerated electrons obeying a power law are injected into the blob uniformly, 
and they are cooled by synchrotron emission and inverse Compton scattering. 



2.1. Kinetic Equations 

The equation describing the time evolution of the electron number spectrum in the blob 
is given by 
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where 7 is the Lorentz factor of electrons and n e {^) is the number density of electrons per 
unit 7. We assume that nonthermal electrons with a power law are injected at a rate: 

9(7) = go7~ P exp(-7/7 max ) for 7 > 7 min , (2) 

where go is the normalization. The injection is assumed to continue during a simulation with 
the rate given above. The energy-loss rate of electrons is denoted by (d^ / dt)\ oss , which is 
due to synchrotron radiation and inverse Compton scattering. The synchrotron-loss rate is 
calculated by the formulation given by Robinson & Melrose (1984) for mildly relativistic 
electrons and by Crusius & Schlickeiser (1986) for relativistic electrons. The cooling rate 
of inverse Compton scattering is calculated by Coppi & Blandford (1990); although the 
EC process is anisotropic, we use the formulation for isotropic scattering. The value of the 
electron-escape timescale, t ej esc, is set to be 3R/c, assuming that the escape of electrons from 
the jet is likely due to advection. We use a method developed by Xiao, Yabe, & Ito (1996) 
for numerical calculations. 

To obtain the isotropic photon field in the blob frame, i.e., except the EC component, 
we solve the following equation: 

dn ph (e) _ n ph (e) 
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where n p h(e) is the photon number density per unit energy e; Compton scattering and 
synchrotron emission are denoted by nc(s) and n cm (e), respectively; n a bs(^) is absorption due 
to the self-absorption of synchrotron process and pair production; and n p ] 1 (e)/t ltesc denotes 
the escape of photons (we set t 7)CSC = R/c). The EC component is separately calculated 
following Georganopoulos et al. (2001; see also Dermer & Schlickeiser 2002) because the 
external photons are highly anisotropic in the blob frame. We assume that in the rest frame 
of the central black hole, the external radiation is isotropic with energy density u cxt and 
monochromatic with energy e cxt . In the blob frame, the comoving energy density of the soft 
photons is u' cxt = w cxt r 2 (l + where = (1 - r~ 2 ) 1/2 (Dermer & Schlickeiser 2002). 

On the other hand, the energy of the external photons is e' ext Te cxt . 

The observed radiation is boosted by the Doppler effect characterized by the beaming 
factor V = [r(l — Pr cos o bs)] _1 , where 9 ohs is the angle between the jet propagation and the 
line of sight (Blandford & Konigl 1979). In this Letter, we assume V — V. The parameters 
to be determined are T, 7 min , 7 max , p, q Q , B, R, e cxt , and u ext . 

The comoving quantities are transformed back into the observer's frame depending on 
the beaming factor and the redshift. The observed photon energy e obs = sV/(l + z), and 
the time duration dt h s = dt (1 + z)/T>. 

3. Numerical Model of a Steady State 

We use the observational data in Hartman et al. (2001). In particular, the spectral 
energy distribution (SED) observed between 1996 January 16 and 30 (P5a in their paper) is 
used for our model of a steady state. 3C 279 is in a preflare stage in this period and is not 
necessarily in a quiescent sate. However, we fitted the data assuming that the spectrum is 
approximated by a steady state solution of the kinetic equations. 

In Figure 1, our model is plotted with the observed data. We obtain the following values 
for the parameters: r = 25, B = 0.3 G, R = 7 x 10 16 cm, p = 1.8, 7 min = 4, 7 max = 1.8 x 10 3 , 
the injection rate 2.3 x 10 _4 cm _3 s _1 , e cxt = 50 eV, and u cxt = 2 x 10" 4 ergs cm" 3 . The 
spectrum below ~ 10 15 Hz is produced by synchrotron radiation, although the radio emission 
is probably from the outer regions of the jet. The radiation between 10 16 and 10 19 Hz is 
from synchrotron-self-Compton scattering, and the 7-rays above 10 20 Hz are from the inverse 
Compton scattering of external soft photons. As shown in Figure 1 by the dashed line, the 
photons scattered by electrons with ~ 7 min appear around 10 20 Hz. Thus the flux around 
10 20 Hz is sensitive to the value of 7 min . However, because the data in this regime are mostly 
upper limits, the constraint on 7 min is weak. 
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Based on our numerical calculations, it is found that mmh/mb = 36, where Wkin is the 
energy density of nonthermal electrons and «g is the magnetic energy density We also find 
that the Poynting power is given by Lp oy = 7tR 2 cubT> 2 = 1.03 x 10 45 ergs s _1 and that the 
kinetic power is given by L^ iQ = TiR 2 cu\a U 'D 2 = 3.69 x 10 46 ergs s -1 . The dominance of the 
kinetic energy of nonthermal electrons over that of magnetic fields is the same as that for 
TeV blazars such as Mrk 421 and Mrk 501 (KTK). 

In Figure 2, the electron spectrum for the SED shown in Figure 1 is presented. It is 
found that the electron spectrum is rather flat in the 7~7 2 n e (7) plot. This is because the 
cooling of high-energy electrons is inefficient, owing to the K-N effects (Blumenthal 1971). 
Note that the external photon energy in the blob frame is e' ~ Te cxt and that the K-N 
effects become effective for electrons with 7k-n ~ m e c 2 /e' ~ 400, which is less than 7 max - 
The spectral shape is contrary to the models often employed to fit the observed data of 
blazars, in which n e {^) oc 7" p below a break energy 7b r and n e (^) oc above 7b r ; the 

value of 7br is determined by t coo i = £ e ,esc, where t coo i is the electron cooling time. Nominally, 
7br is about 30. As seen in Figure 2, the actual spectral shape is much different from this 
conventional expectation. The flat 7-ray spectrum in 10 21 — 10 24 Hz in the v-F v plot is the 
result of these K-N effects. 

The optical depth for photon absorption against electron-positron pair production is 
calculated (Gould & Schreder 1967). The optical depth is much smaller than unity and 
7-rays escape without absorption. 



4. Flare 

The time variability of emission spectra from blazars is quite common. 3C 279 also 
shows time variation from radio to 7-rays (e.g., Wehrle et al. 1998). To understand the jet 
formation mechanisms and particle acceleration processes, it is important to understand the 
mechanism that causes the time variability. Sikora et al. (2001) and Spada et al. (2001) 
show how flare light curves behave according to the shock-in-jet models, which assume that 
blobs with different speeds collide and that the kinetic energy of the bulk motion is dissipated. 
Ballo et al. (2002) interpreted the time variability of 3C 279 as being the result of various 
values of T, although they assumed a broken power-law spectrum of electrons without solving 
the kinetic equation. 

We apply the internal shock model, by which Takahara et al. (2003) will explain the 
time variability of Mrk 421, to a flare observed from 3C 279. Without considering the details 
of the collision process, we simply assume that a coalesced blob has a bulk Lorentz factor Tf 
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and that various physical quantities scale with T. We examine whether or not the flares are 
explained by a change in T, keeping a constant opening angle for the jet cone. The distance 
of the collision location from the base of the jet, d, scales proportional to T 2 , while the lateral 
size of the observed region behaves as R oc T, which is roughly equal to the shell thickness 
in the comoving frame, because the angle from the line of sight is limited by ~ The 
soft photon energy density changes as u ext oc T 2 d~ 2 oc T~ 2 . We assume that the kinetic and 
Poynting powers scale proportional to T 2 so that the total number flux remains constant. 
Thus, the injection rate in the comoving frame, q, scales proportional to rf~ 2 i? _1 oc T~ 5 , and 
the magnetic field scales proportional to T~ 2 . We also assume 7 min and 7 max do not depend 
on T. The value of e' cxt scales proportional to T because of the Doppler effects. Using this 
scaling, we calculate a flare spectrum. Assuming that the spectrum shown in Figure 1 is 
in a quiescent state, the spectrum observed between 1996 January 30 and February 6 is 
fitted. Here the bulk Lorentz factor in Figure 1 is denoted as T q , and that of a shocked 
blob is given by Tf. When Tf = 1.6T q = 40, the 7-ray spectrum is fitted fairly well by our 
model (Figure 3), although there are some minor discrepancies: The slope in the optical- UV 
regime is slightly different; the flux between 10 19 and 10 20 Hz is underestimated; and the 
flux between 10 20 and 10 21 Hz is overestimated. 



5. Summary and Discussion 

We used the kinetic equations of photons and electrons to calculate emission spectra 
from 3C 279, assuming that the emission is from a relativistically moving blob almost along 
the line of sight. Our numerical solution shows that the 7-rays are produced by inverse 
Compton scattering of external soft photons, which is consistent with previous work (e.g., 
Inoue & Takahara 1996; Hartman et al. 2001; Ballo et al. 2002). However, we found 
that the electron spectrum is not a broken power law with a steeper spectrum above a 
break energy. Although electrons are cooled before escape (except those in the lowest energy 
range), the cooling efficiency for high-energy electrons decreases because of the K-N effects. 
As a result, a flat GeV 7-ray spectrum is obtained in the v-vF v plot. It should be noted that 
there is an alternative model of GeV emission that assumes different soft photon sources, i.e., 
the accretion disk and broad line regions (Hartman et al. 2001). In their model, because 
of the different temperatures of soft photons, a broad GeV emission is formed by inverse 
Compton scattering. 

Our result of £ ex t = 50 eV seems to be a little high compared with the conventional value 
of around 10 eV. However, the spectra of ionizing radiation that form the broad line clouds 
are not directly constrained observationally but should extend to energy high enough to 
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multiply ionize atoms of heavy elements. From a theoretical point of view, even the emission 
from the accretion disk is known to deviate from the blackbody because of scattering effects. 
Thus, our choice may be suitable. Recently, it was suggested that IR radiation from dust 
might be important, if the energy dissipation occurs mainly far away from the central region 
(Blazejowski et al. 2000; Sikora et al. 2002); although we have not examined this case, the 
K-N effects will be weaker for this component. 

The kinetic energy density of nonthermal electrons in the blob is an order of magnitude 
larger than the magnetic energy density, which is the same result as for TeV blazars (KTK). 
As a result, the particles transfer more energy in the jet from the central region to the outer 
region than the magnetic fields do. 

We demonstrated that the flare observed between 1996 January 30 and February 6 is 
explained fairly well by the internal shock model with simple scaling laws described by a 
change in the bulk Lorentz factor. It is remarkable that the flare is fitted with those simple 
assumptions. 

It was pointed out by Sikora & Madejski (2000) that the bulk motion of cold elec- 
trons/positrons in the jet scatters the external photons, resulting in observed emission peak- 
ing at energy £bc ~ r 2 e cxt . Our model with T = 25 and e cxt = 50 eV implies 5bc ~ 31 
keV (z/bc ~ 7.5 x 10 18 Hz). The observed luminosity of bulk Compton (BC) scattering is 
estimated as 



™* r 7(r^ r2 h dl/ 



» 3.0 x 10» (1,) 2xlo -!:; gBcm - 3 ( 7xl oi« cm ) "cergss"', 

where o"t is the Thomson cross section, n c is the number density of cold electrons/positrons, 
and the integration is done over the blob volume. For 3C 279 with luminosity distance 
di = 2.4 x 10 3 Mpc, where the Hubble constant H = 75kms~ 1 Mpc~ 1 is assumed, vF v at 
£bc is ~ -^bc/(4vtc?|) ~ 7 x 10~ 14 n c /n e ergs cm _2 s _1 , where n e = 1.65 x 10 4 cm~ 3 is the 
number density of nonthermal electrons shown in Figure 2. The value of n c /n e is unknown, 
and unless it is about 100 or more, the BC spectrum is not observable; note that the observed 
value is ~ 10~ n ergscm -2 s _1 at £bc ~ 31 keV. Regions between the base of the jet and the 
blob may contribute to BC scattering because more cold electrons may exist there than 
in the blob. Let us assume that the bulk Lorentz factor is constant and that the number 
density of the cold electrons behaves as n c oc £~ 2 for £ > £ , where £ is the distance from the 
central region and £q is the critical distance where the bulk Lorentz factor saturates or the 
jet becomes optically thin. If we assume that u ext scales proportional to £~ 2 , the predicted 
value of L B c is d/£ times larger than that of equation (4), noting that n c dV is measured 
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in the comoving frame of the blob. Even if there are no cold electrons is the blob, in the 
inner-jet regions there should exist a corresponding number of cold electrons to relativistic 
electrons in the blob. Thus, we obtain the constraint of > 0.01c?. If n c /n e in the blob is 
more than 100, too much BC emission is predicted, which probably means that u cxt does not 
increase much with a decrease in £. In this case some modifications of our model are needed. 

This work has been partially supported by Scientific Research Grants (M.K. and F.T.: 
13440061; M.K.:15037210; F.T.: 14079205 and 14340066) from the Ministry of Education, 
Culture, Sports, Science and Technology of Japan. 
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Fig. 1. — SED observed between 1996 January 16 and 30 with our model (solid line): The 
data are from Hartman et al. (2001). The values of the parameters are shown in the figure. 
The spectral component of EC scattering is shown by the dashed line. 
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Fig. 2. — Electron spectrum for the SED shown in Figure 1. The vertical dashed line shows 
7 min . The spectrum is flatter above 7 ~ 100 because of the K-N effects. 
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Fig. 3. — SED of the flare between 1996 January 30 and February 6. The solid line is 
calculated assuming Yf/Y q = 1.6. The spectral component of EC scattering is shown by the 
dashed line. 



